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SURFACE  RUNOFF  SUPPLY  ESTIMATES  BASED  ON  SOIL-WATER 
MOVEMENTS  AND  PRECIPITATION  PATTERNS 

By  Leonard  Schiff,  hydraulic  engineer,  Division  of  Drainage  and  Water  Control,  Soil 

Conservation  Service — Research 


THE  PROBLEM  AND  OBJECTIVES  OF  THE  STUDY 

The  infiltration  method  of  estimating  runoff  utilizes  infiltration  curves  developed 
from  hvdrograph  analysis  and  data  obtained  with  inf iltrometers.  The  basic  theory  of  this 
method  has  been  presented  by  R.  E.  Horton,  L.  K,  Sherman,  W.  W.  Horner,  H.  L.  Cook,  and 
others.  Horner  and  Jens  (1)^  have  utilized  infiltration  capacities  and  precipitation  pat- 
terns for  the  development  of  surface  runoff.  More  recently  information  has  been  de- 
veloped on  the  transmission  rate,  (linear  velocity  of  water),  and  percolation  rate,  (rate 
at  which  water  passes  a  specific  soil  level),  within  the  soil  profile.  Considerable  data 
of  this  type  are  now  being  accumulated  under  the  research  program  of  the  United  States 
Deoartment  of  Agriculture. 

The  objective  of  this  paper  is  the  development  of  a  method  of  utilizing  trans- 
mission rates,  antecedent  soil  moisture,  changes  in  soil  moisture,  and  rainfall  patterns 
for  estimating  surface  runoff.     The  method  presented  involves  the  following  concepts: 

1.  Transmission  rates  of  specific  magnitudes  may  be  applied  to  soils  of  reasonably 
similar  hydrologic  characteristics. 

2-  Transmission  rates  may  be  converted  into  corresponding  infiltration  rates. 

3-  An   'jn  place'    supply  of  water   to  the  soil  surface  may  be  estimated.  This 

in  place'  supply  may  be  trans  I  at ed  into  the  hydrograph  of  stream  flow  for 
a  given  sub-basin.  The  hydrographs  for  the  sub-basins  may  be  combined  for  the 
i:ntire  basin.  Est imates  of  the  actual  flow  at  any  po'nt  must  consider  sub- 
surface flow  where  such  flow  is  appreciable. 

EXPERIMENTAL  PROCEDURE 

Rainfall,  runoff,  and  soil  moisture  data  were  obtained  from  two  watersheds  of 
about  1.5  seres,  for  the  period  1946  to  1948,  inclusive.  For  this  period  both  water- 
sheds followed  a  rotation  of  wheat,  meadow,  and  meadow.  Stich  cover  reasonably  pro- 
tected the  soil  surface  against  sealing.  The  soil  on  one  of  these  watersheds  is  Mus- 
kingum silt  loam,  which  is  classified  as  we  11 -drained.  The  soil  on  the  other  water- 
shed is  Keene  silt  loam,  which  is  classified  ss  slowly  permeable. 

Using  infiltration  curves  to  represent  the  supply  of  water  to  the  soil  and  changes 
in  ^ oi 1 -moi s ture  content,  transmission  and  percolation  rates  were  determined  for  Mus- 
kingum silt  loam  and  Keene  silt  loam.  Methods  and  results  of  analyses  were  presented 
by  Schiff  and  Dreibelbis  (2,  3). 

Transmission  and  percolation  rates  were  also  obtained  by  testing  a  few  hundred  core 
samples  of  soil  from  Muskingum  silt  loam  and  Keene  silt  loam.  These  soil  cores  were 
3x3  inches.  Water  was  applied  in  drop  form  close  to  the  surface  of  the  soil,  allowing 
the  drops  to  fall  at  different  points.  Drops  were  applied  rapidly  enough  to  insure 
maximum  rates  of  water  movement  through  the  soil  cores.  The  tim'i  required  for  water  to 
move  through  the  soil  cores  and  appear  at    the  bottom  was  determined.     Velocities  or 


^Italic  numbers  in  parentheses  refer  to  Literature  Cited,  p.  18- 
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t  ransrr.  j  5:  s  i  on  ratp<5  wer^'  then  comouteH.  The  apnJication  of  water  was  continupr*  arr' 
amoMnts  of  water  applied  to  the  surface  of  the  cores  and  amovints  passinP  throuph  the 
cores  were  determined  from  time  to  tim..  The  percolation  rate  was  considered  to  he  the 
rate  of  water  apolication  when  the  volume  of  percolate  equaled  the  volun;e  of  water  ap- 
plied to  the  surface  over  a  cjiven  oeriod  of  time.  This  occiirred  when  the  soil  was  at 
or  near  saturation  The  field  soil-moisture  content  was  also,  determined  for  ench  core. 
*Vthods  and  results  of  analyses  were  presented  hy  Schiff  and  Dreihelhis  (4). 

Terms  and  symbols  used  are  listed  belo"'; 

P  --Accumulated  ground  rein  fait 
i  .-  -Intensi  ty  of  ground  rainfall. 

—A 

Infi  I  trat  ion-  The  passafie  of  water   throufih   the  surface  of  the  soil   into  the 
soi  I  mass. 


/p-  -  In  fi  I  trat  inn  capacity  or  rtotential  infiltration  rate,  that  is,  the  rate 
at  which  infiltration  would  take  place  at  any  instant  were  the  sunply 
to  equal  or  exceed  this  capacity. 


Q^-  -Accumulated  surface  runoff. 

TjQg  and  T^,-  -Theoret ical  accumulated  surface  water  supply  to  runoff,  that 
i  s ,  the  supply  equal  to  'a'^p  f'"'"'  subsur  face  f  low  contr  ibut  ions . 
(Tj  and  refer  to  two  different  methods.) 

_S--/'ore  space  or  -  toraie  space  available  in  the  soil  at  a  /Jfven  time. 

Transmission-  The  movement  of  water  downward  through  a  given  depth  within  the 
soil  mass,  that  is,  a  linear  velocity  of  water  through  a  given  depth  of 
so  I  / . 


^--Transmission  rate,      t^^     is   the   transmi  ss  ion  rate  in   the     topsoil , 
t is  the  transmission  rate  in  the  subsoil. 

Percolation-  The  water  passing  a  given  level  within  the  soil  mass,  that  is, 
a  depth  of  water  over  a  given  area  passing  a  specific  soil  level  in  a 
specific  period  of  time.  It  is  equal  to  the  area  of  pore  space  multi- 
plied hy  the  transmission  rate  through  this  pore  space  for  a  given 
period  of  time,  St. 

^-Percolation  rate. 

STEPS  IN  ESTIMATING  SURFACE  RUNOFF 


Steps   in  estimating  surface  runoff  are 

1.     Prepare  a  table  showing  the  theoretical 
various  soi l-moisture  contents  to  absorb 


outlined  below  and  detailed  later: 

transmission  rates  required  in  soil  at 
various  infiltration  rates. 


2.  In  this  table  place  a  line  under  the  transmission  rate  nf  the  topsoil,  a  tine 
under  the  transmi  ssion  rate  of  the  subsoil,  and  other  portions  nf  the  soil  pro- 
file if  necessary ,    for  the  particular  soil  studied. 


3. 


Using  this  table  prepare  potential^  infi  Itration  curves,  for  selected  ante 

cedent  soi I -moisture  contents  and  changes  in  moisture  contents . 
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4.  These  curves  may  be  applied  to  an  ant  icipated  storm,  a  par  ticu  lai  storm, 
or  to  a  specific  storm  pattern  that  may  be  expected  in  any  locali  ty.  The  dif- 
ference between  the  /jround  intensity  bargraph  and  the  potenti  al  infiltration 
curve,     'g'^p'  place'   surface  supply  of  water  to  runoff. 

5.  This  'in  place'  supply  of  water  may  be  translated  into  the  surface  runoff  hydro- 
graph  of  stream  flow  for  a  given  sub-  basin. 

6.  The  surface  runoff  hydrographs  for  the  sub-  basins  may  be  combined  for  the  entire 
basin  and  when  further  combined  with  the  s'ubsur  face  flow  hydrograph  may  be 
t ranslat ed  into  stream  flow. 

Steps  5  and  6  are  not  di«!cussed  herein.  They  involve  the  translation  of  ^i'^p  i"' 
to  stream  flow,  such  as  through  surface  detention  and  channel  storage  coirputat ions . 

To  avoid  numerous  repetition  of  units  in  the  following  discussion,  soil -moi sture 
content  is  expressed  in  inches  of  water  per  inch  of  soil.  Field  capacity  is  nearly 
enough,  0.30  inch  of  water  per  inch  of  soil  for  the  topsoil  and  subsoil  of  Muskinguir  silt 
loam  and  the  topsoil  of  Keene  silt  loam  and  0.36  inch  for  the  subsoil  of  Keene  silt  loam 
Saturation  is  nearly  enough,  0.50  inch  for  the  topsoil  and  subsoil  of  both  Muskingum  and 
Keene  silt  loam.  The  topsoil  is  the  deoth  of  soil  from  0-7  inches  and  the  subsoil  is 
the  depth  of  soil  from  7-14  inches. 

Step  1  -  Studies  (2),  (3),  (4)  and  by  others,  show  that  a  considerable  rate  of 
ground  rainfall  establishes  somewhat  of  a  'saturated  front'  that  moves  downward 
through  soils  of  the  type  considered  herein.  The  infiltration  rate  is  a  func- 
tion of  the  transmission  rate  within  the  soil.  Table  1,  page  4,  presents 
theoretical  transmission  rates  required  in  soil  at  a  given  soil-moisture  con- 
tent to  absorb  various  infiltration  rates.  For  example,  the  transmission  rate 
of  20  inches  per  hour,  shown  under  the  infiltration  rate  of  4.0  inches  per  hour 
and  opposite  a  soil-moisture  content  of  a  0.30  or  0.20  inch  of  available  stor- 
age space  per  inch  of  soil  was  com.puted  as  follows:  It  requires  3  minutes  to 
fill  0.20  inch  of  storage  space  with  water  at  an  infiltration  rate  of  4.0 
inches.  Since  0.20  represents  the  available  storage  space  in  an  inch  of  soil, 
water  must  have  moved  1  inch  in  3  minutes  and  the  transmission  rate  is  20 
inches  per  hour.  In  other  words,  0.20  inch  of  storage  space  multiplied  by  a 
transmission  rate  of  20  inches  per  hour  equals  an  infiltration  rate  of  4.0 
inches  per  hour.  Other  values  in  table  1  were  computed  accordingly.  If  water 
did  not  move  as  a  saturated  front  it  would  have  moved  more  than  1  inch  and  the 
transmission  rate  would  have  been  greater.  This  depends  on  the  characteristics 
of  the  soil. 

The  percolation  rate  at  a  given  level  within  the  soil  is  equal  to  the  in- 
filtration rate  if  the  soil  is  saturated  above  that  level.  Prior  to  satura- 
tion the  percolation  rate  is  equal  to  the  infiltration  rate  minus  the  rate  of 
exhaustion  of  pore  space  above  the  level  considered. 

Step  2  -  The  stepped  solid  line  in  table  1  has  been  placed  under  20  inches  per  hour, 
the  maximum  transmission  rate  generally  found  for  Muskingum  silt  loam  and 
Keene  silt  loam  at  Coshocton,  Ohio  (2),  (3),  and  (4).  The  dot  and  dash  line 
represents  Muskingum  silt  loam  subsoil  and  has  been  placed  under  6.7  inches  per 
hour.  The  dotted  line  represents  Keene  silt  loam,  subsoil  and  has  been  placed 
under  5.0  inches  per  hour.  Subsoil  values  taken  from  references  (2)  and  (3) 
were  lower  than  those  obtained  by  ^ore  sampling  (4).  Values  shown  in  reference 
(2)  and  (3)  were  obtained  by  using  infiltration  rates  as  the  supply  of  water  to 
the  soil  and  frequently  this  supply  was  inadequate  for  maximum  transmission 
rates  in  the  subsoil.  Values  showp  in  reference  (4)  are  higher  for  the  subsoil 
and  were  obtained  by  the  application  of  an  adequate  supply  of  water  to  soil 
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cores.  Obviously,  the  ability  to  correctly  locate  the  stepped  lines  will  de- 
termine the  accuracy  of  relationships  between  soil  moisture  and  infiltration 
rates.  Differences  in  transmission  ^p.tes  do  mt  appear  Rieat  between  the  sub- 
soils of  Muskingum  silt  loam  and  K -ene  silt  loam.  Ho\?'ever,  the  differences 
that  do  exist  as  weJl  as  the  greater  swelling  ot  colloids  in  the  subsoil  of 
Keene  silt  loam,  seem  to  cause  considerable  differences  in  hydrologic  per- 
formance. This  is  evidenced  by  the  higher  mois^^ure  content  of  Keene  silt  loam 
antecedent  to  a  storm,   a  major  reason  for  greater  runoff  from  Keene  silt  loam. 

Step  3  •  For  Muskingum  silt  loam  and  Keene  silt  loam  the  water  movement  within  the 
soil  profile  affecting  the  supply  of  water  available  for  surface  runoff  takes 
place  within  the  topsoil  and  the  subsoil  (2),  (3).  When  the  soil  is  below 
field  capacity  and  the  supply  of  water  is  sufficient  for  infiltration  to  occur 
near  capacity  rates,  water  moves  rapidly  downward  until  it  reaches  the  sub- 
soil. Under  such  conditions  when  water  reaches  the  subsoil  the  soil -moisture 
content  of  the  topsoil  will  be  about  0.46.  When  the  topsoil  is  above  field 
capacity  water  will  move  into  th^  subsoil  regardless  of  an  approaching  wet 
front.  HoT"ever,  for  most  practice  1  problems  of  surface  runoff  it  may  be  as- 
sumed that  water  wi  1  enter  the  subsoil  when  infiltration  exists  near  capacity 
rates  and  the  soil -moisture  content  of  the  topsoil  is  equal  to  or  greater  than 
0.46.  Tn  view  of  these  consideration-  the  soil -moisture  contents  of  the  top- 
soil  and  subsoil  overlap  at  0.46  on  the  horizontal  scale  in  figure  1,  page  6. 
When  the  soil -moisture  content  of  the  topsoil  exceeds  0.46  the  low  transmission 
rates  of  the  subsoil  become  the  controlling  factor  and  surface  runoff  becomes 
significant.  Since  one  unit  on  the  horizontal  scale  represents  an  increase  in 
moisture  content  of  0.01  inch  of  water  per  inch  of  soil  or  0.07  inch  pain  in 
the  topsoil,  two  units  are  necessary  for  both  the  topsoil  and  subsoil  to  gain 
0.01  inch  of  water  per  inch  of  soil  or  0.14  inch  gain  in  both  the  topsoil  and 
subsoil.  The  scale  was  adiusted  accordingly  in  the  overlapping  section.  For 
this  overlapping  section  the  division  of  the  total  gain  would  be  more  accurate 
if  based  on  the  ratio  of  the  storage  space  depletion  rate  of  the  topsoil  to  the 
subsoil.     Such  refinements  are  unnecessary  for  most  practical  solutions. 

Table  1,  page  4,  was  used  to  prepare  the  potential  infiltration  rate  curves 
for  various  antecedent  soil-moisture  contents  shown  in  figure  1.  Plot- 
ting potential  infiltration  ra*es  against  soi  i -moi sture  results  in  straight 
lines  of  different  slope,  one  for  the  topsoil,  one  while  the  remaining  storage 
space  is  being  exhausted  in  the  topsoil  as  water  moves  through  the  subsoil,  and 
one  for  the  subsoil  after  the  topsoil  is  saturated.  However,  in  figure  1  the 
potential  infiltration  rates  for  different  soil-moisture  contents  were  plotted 
an  smooth  curves,  using  the  straight  lines  as  guides. 

As  the  soil-moisture  content  increases,  the  effective  cross- sectional 
area  through  which  water  flows  downward  decreases,  although  transmission  rates 
or  velocities  ma-'  remain  approximatejy  uni 'orm  and  fairly  high.  Diminishing 
effective  areas  of  cross  section  multiplied  by  the  fairly  constant  velocity  of 
flow  resutt  in  diminishing  rate."  of  flow  through  the  soil  and  consequently 
diminishing  infiltration  rates.  After  water  begins  to  move  slowly  into  a  sub- 
soiJ  of  lower  transmission  rates,  the  remaining  storage  space  fills  in  the  top- 
soil  and  the  topsoil  approaches  saturation.  The  lower  rates  of  exhaustion  of 
this  remaining,  and  generally  smaJl  volume  of  storage  space  in  the  topsoil, 
and  lower  transmission  rates  in  the  subsoil  usually  are  reflected  by  the  sec- 
tion of  greatest  curvature  in  the  infiltration  capaci ty- rate  curve  when  plot- 
ted against  time.  Values  for  this  section  were  obtained  frow  the  subsoil  steps 
in  table  1.  The  infiJtration  curve,  after  the  topsoil  is  saturated,  is  mainly 
a  function  of  water-movement  rates  in  the  subsoil.  The  infiltration  rate  is 
further  reduced  as  storage  space  is  exhausted  in  the  subsoil  and  probably  ap- 
proaches a  minimum  as  the  swelling  of  colloids  beccmes  complete. 


ANTECEDENT  SOIL  MCTSTURE 
TOPSOIL  AND  SUBSOIL 
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Sections  makint;  up  the  Ruidc  lines, .from  which  the  potentiaJ  infiltration 
rate  curves  are  drsv/n,  would  not  be  straight  lines  if  transmission  rates  varied 
with  progressive  changes  in  soil  moisture.  However,  transmission  rates  couJd 
be  constant  for  a  given  horizon  in  the  soil  profile  and  storage  space  exhausted 
in  such  a  fashion  as  to  result  in  the  curvature  of  various  sections  of  the  in- 
filtration rate  curve. 

The  curves  shown  in  fig-jre  1,  page  6,  were  prepared  for  the  soiJs  at  field 
capacity,  0.40,  0.45,  and  for  combi  n.at  i  ons  representing  antecedent  moisture 
conditions  existing  for  the  storms  of  Tune  16,  1946,  and  April  14,  1948.  Such 
combinations  are  marked  on  each  curve.  For  example  the  top  curve  is  marked 
MSI.  .21,  .24,  which  represents  Musking>jm  silt  loam  at  a  oil-moisture  content 
of  0.21  for  the  topscil  and  0.24  for  the  subsoil.  This  was  the  soi I -moi sture 
condition  for  the  storm  of  June  16,  1946,  on  watershed  109.  Numerous  curves 
may  be  prepared  in  advance  covering  a  number  of  likely  soil-moisture  condi- 
tions. 

Step  4  -  This  step  determines  the  supply  of  water  to  the  soil  surface  using  the 
potential  infiltration  rate  curves  shown  in  figure  1  and  the  rainfall  pattern. 
This  supply  of  water  will  be  equal  to  surface  runoff  if  depression  storage  is 
negligible.  Information  from  the  previous  steps  may  be  developed  prior  to 
storms  to  which  it  is  to  be  applied.  Two  met!-.ods  have  been  developed  to  apply 
this  information.  The  first  method,  Tj,  is  more  rapid  but  not  as  accurate 
as  the  second  method,  Method    Tj     assumes  that  all  of  the  ground  rainfall 

enters  the  soil,  is  uniformly  distributed  in  the  topnoil  until  the  moisture 
content  of  the  topsoil  is  0.46  inch,  is  then  equally  divided  between  the  top- 
soiJ  and  subsoil  until  the  topsoil  is  saturated,  and  is  then  uniformly  dis- 
tributed in  the  subsoil.  When  the  subsoil  becomes  saturated  the  minimum  'infil- 
tration rate  occurs.  Assuming  that  all  of  the  ground  rainfall  enters  the  soil, 
and  omitting  any  correction  for  surface  runoff,  obviously  assumes  that  more 
water  enters  the  soil  than  is  actually  the  case.  However,  this  is  somewhat 
compensated  for  by  assuming  uniform  distribution  of  water  in  the  subsoil  after 
the  topsoil  is  saturated.  Actually,  water  moves  slowly  through  the  subsoil, 
and  the  upper  portion  of  the  subsoil  is  wetter  than  the  lower  portion  for  a 
considerable  length  of  time. 

Method    Tj    is  illustrated  in  figure  2,  page  8,   for  the  storm  of  June  16, 
1946,  on  watershed  123,  Keene  silt  loam,  slowly  permeable.     Plot  the  accumu- 
lated ground  rainfall,     '*g>     and  the  intensity  of  ground  rainfall  or  intensi- 
ty bargraph,     ig.     The  intensities  of  ground  rainfall  may  be  placed  along  the 
curve  corresponding  to  the  slopes  of  the  elements  of  the  curve,  elim* 

iiiating  the  need  for  the  intensity  bargraph.  Next,  the  soil -mdvsture  ordinates 
for  the  topsoil  and  subsoil  are  established.  The  antecedent  soil-moisture 
content  for  the  storm  of  June  16  was  0.29  for  both  the  topsoil  and  subsoil. 
The  value  of  0.29  for  the  topsoil  is  placed  opposite  the  ground  rainfall  of 
0.00  inch  on  the  vertical  sc.ile.  For  each  0.07  inch  of  ground  rainfall  the 
soil-moisture  content  of  each  inch  of  topsoil  is  increased  by  0.01  until  the 
value  of  0.46  is  reached,  in  this  case  opposite  1.19  inches  of  ground  rainfall. 
The  antecedent  value  of  0.29  for  the  subsoil  is  also  placed  opposite  1.19 
inches,  since  at  this  point  water  begins  to  enter  the  subsoil.  For  each  ad- 
ditional 0.14  inch  of  ground  rainfall,  the  soil -moisture  content  of  each  inch 
of  topsoil  is  increased  by  0.01  until  the  topsoil  is  saturated,  in  this  case 
opposite  1.75  inches  of  ground  rainfall.  After  this  time,  for  each  additional 
0.07  inch  of  ground  rainfall,  the  soil -moisture  content  of  each  inch  of  sub- 
soil is  increased  by  0.01  until  the  subsoil  is  saturated  and  the  minimum  infil- 
tration rate  established. 
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Plot  potential  infiltration  rates  on  figure  2,  page  8,  obtained  from  curve 
KSL,  .29,  .29  in  figure  1,  page  6,  corresponding  to  the  ordinate  of  soil-mois- 
ture values  established  on  figure  2  and  the  accumulated  ground  rainfall,  Pa 
For  example,  for  a  soil-moisture  content  of  0.35  for  the  topsoil,  curve  KSL  .29, 
.29  of  figure  1  shows  a  potential  infiltration  rate,  fp,  of  2.9  inches  per 
hour,   shown  as     ^l^p  method    Tj     on  figure  2.     The  value  2.9  is  plotted 

above  the  intersection  of  the  soil-moisture  ordinate  of  0.35  and  the  Pg  curve. 
Similarly,  plot  sufficient  other  fp  points  corresponding  to  other  values  of 
soil— moisture  content  to  establish  the  curve  Tjfp. 

Intensities  of  ground  rainfall,  ig,  in  excess  of  the  ^f'p  curve,  yield 
water-supply  increments  available  for  surface  runoff.  For  example,  from  9:17 
p.  m.  to  9:21  p.  m.  ig  is  4.05  inches  per  hour  and  the  average  ^l^p 
this  period  of  time  is  1.54  inches  per  hour,  a  difference  of  2.51  inches  per 
hour.  In  this  period  of  4  minutes,  a  rate  of  2.51  inches  per  hour  is  equal  to 
0.167  inch  of  water  supply.  Starting  where  the  fp  curve  first  intersects  the 
ig  bargraph  compute  each  increment  of  water  supply  and  accumulate  the  incre- 
ments to  establish  TjQg,  the  theoretical  accumulated  water  supply.  If  there 
is  little  depression  storage    TjQg    may  be  considered  surface  runoff. 

At  midnight  of  June  16,     ^iQg  1-32  inches  as  compared  to  1.15  inches 

of  actual  surface  runoff,  Q^.  It  is  believed  that  some  of  this  difference  is 
due  to  the  assumption  that  all  of  the  ground  rainfall  entered  the  soil.  Other 
differences  are  interesting.  For  example,  since  the  TjOg  curve  represents 
the  total  available  water  supply  on  the  watershed  at  a  given  time  for  a  given 
increment  of  rainfall  it  starts  earlier  and  steeper  than  the  Qg  curve. 
Initial  runoff  shown  by  a  runoff-measuring  flume  usually  comes  from  the  area 
in  the  immediate  vicinity  above  the  flume.  Obviously,  the  lag  prevails 
throughout  the  runoff  period  as  does  the  ironing  out  of  the  increments  of  water 
supply.  These  differences  reflect  the  influence  of  the  physical  factors  of  the 
surface  of  the  watershed  upon  the  movement  of  the  water  supply  to  the  measuring 
flume,  where  it  is  measured  as  surface  runoff. 

Method  T2  is  also  illustrated  in  figure  2  for  the  storm  of  June  16, 
1946,  on  watershed  123.  This  m.ethod  is  similar  to  method  .  except 
that  it  accounts  for  the  water  entering  the  soil  by  subtracting  the  total 
increments  of  estimated  surface  runoff  at  any  given  time  from  the  accumulated 
ground  rainfall  up  to  the  same  time.  Also,  the  distribution  of  water  within 
the  topsoil  and  subsoil  is  made  more  strictly  in  accordance  with  soil-water 
movements.  Table  2,  page  10,  shows  how  the  accumulated  surface  runoff  curve, 
T^s  °^  figure  2,  was  obtained.  A  few  fp  points  are  determined  to  locate 
and  shape  the     fp    curve,  shown  as    ^2^p  '"^t^'o^    "^2'     Prior  to  its  inter- 

section with  the  intensity  bargraph.  Column  1  of  table  2  shows  such  points 
to  be  taken  at  9:02M,  9:05,  9:07,  and  9:08H  p.  m.  when,  between  each  of  these 
times,  increments  of  ground  rainfall  were  0.07  inch,  column  2.  For  each  of 
these  time  entries  the  gain  in  soil-moisture  content  of  the  topsoil  is  shown  in 
columns  10  and  11.  Values  of  fp  corresponding  to  these  soil -moisture  con- 
tents were  obtained  from  figure  1  and  are  3.11,  2.91,  2.72,  and  2.54  inches 
per  hour  as  shown  in  column  14.-  A  curve,  through  these  points,  was  extended 
until  it  intersected  the  intensity  bargraph  at  9:08H  p.  m.  This  curve  was 
continued  through  the  intensity  bargraph  section  whose  value  is  2.54  inches 
per  hour.  For  the  period  9:08^  to  9:11  p.  m.  the  difference  between  the  ig 
value  of  2.54,  column  3,  and  the  trial  average  fp  value  of  2.44,  column  4, 
is  0.10  inch  per  hour,  column  5.  This  difference,  for  the  time  interval  of 
2hi  minutes,  column  6,  is  equivalent  to  a  trial  increment  of  surface  runoff  of 
0.004  inch,  column  7.  Column  8  is  the  final  iAPs  of  the  previous  line, 
column  18,  plus  the  trial  of  column  7  for  the  interval  being  considered. 

Subtracting  the  trial  APj    of  0.004  inch  plus  the  final  fAPg,     column  18 
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(in  this  case  0.000),  from  a  Pg  of  0.66  yields  0.656  inch  as  the  total  soil- 
moisture  gain  up  to  that  time,  column  9.  This  entry  of  0.656  is  0.096  inch 
greater  than  the   final  minus  "^AQg    of  0.56,   column  19,  previous  line. 

The  value  0.096  for  the  entire  topsoil  represents  a  gain  of  about  0.014  inch 
of  water  per  inch  of  toosoil,  column  10.  This  gain  of  0.014  when  adder)  to  the 
previous  entry  of  0.37  inch  in  column  11  is  the  total  soil-moisture  content  of 
0.384.  The  zero  in  column  12  indicates  no  moisture  gain  in  the  subsoil,  7-8 
inches.  Consequently,  the  total  soil -moisture  content  of  the  subsoil  rpmains 
the  same,   0.29,  column  13.     The  value  obtained  from  figure  1,  page  6,  for 

a  soil— moisture  content  of  0.384  {column  11)  is  2.27  inches  per  hour,  column 
14.  and  is  considered  a  final  position  point  through  which  the  fp  curve 
passes.  The  final  average  fp  for  the  interval  9:08^  to  9:11  p.  m.  is  the 
2.40  of  column  15,  rather  than  the  trial  value  of  2.44  of  column  4.  The  final 
ig  minus  the  average  /_  Is  the  0.14  of  column  16  rather  than  the  0.10  of 
column  5.  The  final  LOg  is  the  0.006  of  column  17  rather  than  the  0.004  of 
column  7.  The  0.006  of  column  17  is  the  final  i.LQg,  which  when  subtracted 
from  the  Pg  of  column  2  is  the  0.654  of  column  19.  This  procedure  was  con- 
tinued step  by  step  for  the  other  intervals  of  the  intensity  bargraph.  For 
greater  accuracy,  a  long  interval  of  constant  intensity  may  be  divided  into 
parts. 

To  illustrate  the  method  when  water  has  reached  the  subsoil  and  is  de- 
pleting storage  space  in  both  the  topsoil  and  subsoil,  consider  the  interval  of 
intensity  from  9:31  to  9:33  p.  m.  The  value  of  1.254  shown  in  column  9  is 
0.04  greater  than  the  value  of  1.214,  previous  line,  column  19.  The  0.04  was 
proportioned  between  the  topsoil  and  subsoil  on  the  basis  of  their  storage 
space  depletion  rates  (2).  This  resulted  in  0.023  inch  gain  for  the  topsoil 
and  0.017  for  the  subsoil.  The  topsoil  gain  was  considered  to  be  0.003  inch 
per  inch  as  shown  in  column  10  and  0.017  inch  for  the  7-8-inch  depth  of  the 
subsoil  as  shown  in  column  12.  Similarly  the  remainder  of  table  2,  page  10, 
was  completed  and  column  18  was  plotted  on  figure  2,  page  8,  as  T^^Ps'  Total 
surface  runoff,     "^^s'  1.09  as  compared  to  an  actual  surface  runoff,  Q^, 

of  1.15  inches.  If  the  jT^Pg  curve  were  absolutely  correct  the  area  between 
the    7*2^3    and  curves,  after  their  intersection,  could  be  considered  the 

amount  of  subsurface  flow  return. 

Where  the  slope  of  the  T2P5  curve  is  constant  for  a  period  of  time  at  least  equal 
to  the  'time  of  concentration'  of  the  watershed  it  may  be  considered  an  approximation  of 
the  rate  of  runoff.  For  example,  the  rate  of  runoff  based  on  the  slope  of  "^^Ps  ''ver 
the  period  of  time  from  9:17  to  9:31  p.  m.  is  2.31  compared  to  an  actual  rate  of  1.91 
inches  per  hour  occurring  at  9:32  p.  m.  Rates  of  runoff  may  be  estimated  for  slopes  for 
periods  less  than  the  time  of  concentration. 

Figure  3,  page  12,  shows  the  two  methods  for  the  same  storm  on  watershed  109,  Mus- 
kingum silt  loam,  well -drained .  Less  runoff  on  watershed  109  than  on  watershed  123  is 
due  more  to  its  lower  antecedent  soil-moisture  condition  prior  to  a  storm  thsn  to  its 
higher  transmission  rates.  There  are  more  soil  particles  of  colloidal  dimension  in 
Keene  silt  loam  and  greater  swelling  probably  accounts  largely  for  the  poorer  drainage 
during  and  between  storm.s  on  watershed  123. 

Figure  4,  page  13,  shows  T^Pg  for  the  storm  of  April  14,  1948,  on  watersheds  123 
and  109.  Of  particular  note  is  the  shape  of  fp,  which  shows  sudden  drops  with  sudden 
increases  in  rainfall  intensity  and  soil-moisture  content.  There  is  close  agreement  be- 
tween theoretical  and  actual  runoff.  If  there  are  long  periods  in  a  storm  during  which 
no  rain  falls  there  will  be  n  decrease  in  the  soil-moisture  content,  whict>  may  be  esti- 
mated from  soil-moisture  depletion  rates  (3).  Storage  space  and  the  potentiaJ  infiltra- 
tion rate  are  thus  increased.  When  rain  begins  to  fall  again  the  infiltration  rate  will 
be  based  on  the  actual  storage  space  at  the  time  and  will  show  the  recovery  that  seems 
somewhat  difficult  to  establish. 
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These  methods  may  be  appJied  to  an  anticipated  storm,  a  particular  storir  or  to  a 
specific  storm  pattern:  Potential  infiltration  curves  may  be  prepared  for  various  ante- 
cedent soil-moisture  conditions  and  superimposed  over  various  storm  classes  and  patterns 
that  might  occur  in  a  given  locality.  Five  classes  and  six  patterns  have  been  developed 
for  rains  at  Coshocton,  Schif f  (5).  The  classes  refer  to  various  combinations  of  rainfall 
intensities  of  certain  magnitudes  and  the  patterns  refer  to  the  sequence  of  these  in- 
tensities. Seasonal  occurrences  of  these  classes  of  storms  in  the  province  represented 
by  Coshocton,  are  generally  as  follows:  Class  T,  winter;  class  IT,  winter  and  early 
spring;  class  III,  fall;  class  TV,  late  spring  and  summer;  and  class  V,  summer.  The  ad- 
vanced and  delayed  patterns,  highest  intensities  occurring  near  the  beginning  and  end  of 
storms,  respectively,  are  used  for  illustration  herein,  as  all  other  patterns  fall  be- 
tween these.  Whenever  possible,  antecedent  soil-moisture  conditions  should  be  obtained  by 
actual  field  determinations  of  the  soil -moisture  content  rather  than  by  estimates  based 
on  rainfall.  Better  relationships  between  rainfall  and  soil-moisture  content  might  be 
established  in  the  light  of  information  on  rainfall,  percolation,  and  evapo-transpi ration. 

Figures  5  and  6,  pages  15  and  16,  show  fp  curves  obtained  from  figure  1,  page  6, 
for  various  antecedent  soil-moisture  conditions,  applied  to  advanced  and  delayed  patterns 
of  five  classes  of  rainfall.  Resulting  runoff  based  on  method  Tj^  is  also  shown  in 
figures  5  and  6 .  For  storm  classes  I,  TI,  and  III,  an  antecedent  soil-m.oi sture  content 
of  0.40  was  used.  For  stcrm  class  IV  and  V  field  capacity  was  used,  with  the  exception  of 
one  soil-moisture  condition  taken  at  0.21,  0.24  for  Muskingum  silt  loam  and  0.29,  0.29  for 
Keene  silt  loam,  soil-moisture  contents  actually  existing  prior  to  the  storm  of  June  16, 
1946.  A  delayed  pattern  storm  produces  more  runoff  than  an  advanced  pattern  storm  as  the 
high  intensities  fall  on  soil  whose  moisture  content  has  increased  under  the  influence  of 
the  lower  intensities  at  the  beginning  of  the  storm.  Greater  runoff  occurs  on  Feene  silt 
loam  for  reasons  given  previously.  There  is  close  agreement  between  the  actual  runoff  for 
the  storm  of  June  16  1946,  a  class  IV,  Pattern  A  type  storm,  and  theoretical  runoff  for 
such  a  type  storm  and  similar  soil-moisture  content  on  watersheds  123  and  109,  as  shown 
in  figure  6. 

SUMMARY 

This  paper  attempts  to  set  up  a  fundamental  attack  that  involves  two  important  ob- 
jectives. The  first  is  concerned  with  estimates  of  surface  supplies  to  runoff  and  is 
dealt  with  directly.  Potential  infiltration  curves  can  be  used  to  estimate  surface  sup- 
plies to  runoff  by  superimposing  them  over  storm  patterns  of  various  types  that  have  oc- 
curred or  are  expected  to  occur .  Infiltration  curves  based  on  soil-water  movements  should 
give  a  better  understanding  of  the  ranges  and  recoveries  of  such  curves  for  different 
soils  and  land  use.  The  second  objective,  possibly  of  even  greater  ultimate  importance, 
although  not  emphasized  in  this  particular  paper,  is  an  approach  which  will  lead  to 
knowledge  of  the  influence  of  soil  properties  and  land  use  upon  the  hydrologic  performance 
of  watersheds.  Soil-water  movements  represent  the  supply  of  water  to  the  soil  profile 
and  are  important  considerations  in  land  use  and  subsur face- type  flows.  The  actual  rate 
of  flow  at  any  given  point  in  the  drainage  basin  is  a  function  of  subsurface  flow  contri- 
butions as  well  as  the  action  of  the  physical  characteristics  of  the  watershed  surface  on 
the  accumulative  surface  runoff  curve. 

Laws  may  be  developed  from  such  an  approach  which  may  be  applicable  anywhere,  al- 
though certain  factors  may  vary  in  magnitude  based  on  the  soil  and  climatic  province  under 
consideration.  This  paper  attempts  to  set  up  relationships  between  soil -water  movem.ents 
and  surface  runoff  that  can  be  applied  to  other  soils.  However,  it  would  be  necessary  to 
determine  certain  hydrologic  characteristics  of  other  soils  as  was  done  for  the  soils 
mentioned  herein.  Since  the  approach  in  this  paoer  is  believed  to  be  a  fundamental  one  it 
does  not  preclude  soils  which  tend  to  seal  or  crack.  Velocity  or  transmission  rates 
through  such  soils  may  also  be  approximated  and  such  rates  may  be  related  to  the  pore 
space  or  crevices  and  cracks  through  which  major  movements  of  water  downward  occur. 
Seasonal  differences  in  the  hydrologic  characteristics  of  soils  may  also  be  accounted  fot 
in  this  approach  as  consideration  may  be  given  to  seasonal  changes  in  soil -water  move- 
ments. 


ELAPSED  TIME  IN  HOURS  ELAPSED  TIME  IN  HOURS 


FIG.  5 -SURFACE  RUNOFF  ESTIMATES  BASED  ON  SOIL-MKTER  MOVEMENTS  AND  PRS3nW)0N  PATTERNS 
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It  would  seem  desirable  to  determine  the  hydrologic  characteristics  of  major  soils 
of  the  country  that  show  important  differences  in  permeability.  These  soils  could  then 
be  grouped  hydrologically.  Soils  for  which  such  determinations  were  not  practical  could 
probably  be  fitted  reasonably  well  into  the  established  groups.  Such  a  grouping  related 
to  a  grouping  cf  climatic  provinces  would  be  a  powei ful  tool  in  evaluating  land-use  prac- 
tices and  in  problems  involving  runoff. 
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